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Abstract 

The perturbation theory with respect to the potential energy of three particles is considered. 
The first-order correction to the continuum wave function of three free particles is derived. It is 
shown that the use of the collective multipole expansion of the free three-body Green function over 
the set of Wigner D-functions can reduce the dimensionality of perturbative matrix elements from 
twelve to six. The explicit expressions for the coefficients of the collective multipole expansion 
of the free Green function are derived. It is found that the S-wave multipole coefficient depends 
only upon three variables instead of six as higher multipoles do. The possible applications of the 
developed theory to the three-body molecular break-up processes are discussed. 
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I. INTRODUCTION 



The study of angular distributions in processes of three-particle fragmentation is an im- 
portant source of information about the dynamics of many physical objects such as atoms 
and molecules. 

Till recently, the experimental technique allowed one to analyze the angular distributions 
only for charged particles resulting from the many-body fragmentation. In such situations, 
the structure of the angular distributions is determined by the Coulomb force. For example, 
this is realized in the process of two-electron single-photon ionization of Helium atom which 
is quite well studied both experimentally and theoretically [lj]. 

During the last decade, the progress of experimental technique has made it possible to 
investigate the angular distributions of neutral fragments arising in the processes of three- 
atomic molecular break-up. Namely, in {2] the dissociative recombination of water ion into 
the neutral atomic fragments has been studied. In a series of experiments the 
angular distributions in the predissociation of the triatomic hydrogen into three hydrogen 
atoms were investigated. Clearly, the analysis of angular distributions in molecular break-up 
processes could gain significant insight into the nature of chemical forces. 

In molecular physics the dynamics of the interacting atoms is usually described within 
the framework of the Born-Oppenheimer approximation. In this approach the interaction 
of ions is represented by the set of potential energy surfaces (PES) which are, in fact, the 
mean field of the electrons. It is important that, generally, PES cannot be decomposed 
into a sum of binary potentials (i.e. pairwise interactions). Rather, it will contain the terms 
which "entangle" the coordinates of all particles thus making the standard methods (such as 
Faddeev approach Q]) inapplicable. Thus, the development of an adequate theory of many- 
particle molecular fragmentation is a complicated task. Therefore, it is highly desirable to 
have simple yet physically meaningful models of the molecular break-up. 

In the break-up into the neutral fragments often the situation is realized when the kinetic 
energy of fragments prevails over their interaction potential. This opens the possibility 
to apply the perturbation theory to the calculation of the wave function of three-particle 
continuum. In the presented paper the lowest-order perturbation theory with respect to the 
potential energy was applied to the quantum three-body system. Even in this simplest case 
the calculation of the matrix elements leads to twelve-dimensional integrals. However, it 
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turns out that the dimensionality of these integrals can be reduced from twelve to six 
employing the technique of collective multipole expansions. This technique was used in 
where such an expansion was derived for the product of two three-dimensional plane waves 
expz(q! ■ i - ! + q 2 ■ r 2 ). Note that the results presented in the paper allows one to simplify 
the calculation of the matrix elements within the perturbative approximation even if the 
potential contains three-body (i.g. non-binary) terms. 

The paper is organized as follows. In Sec. [Ill the general equations of the perturbation 
theory are derived. They are based on the expression for the Green function corresponding to 
a system of three free particles. The calculation of matrix elements involving the perturbative 
wave functions is considered in Sec. IIHI In that section it is also demonstrated how the 
calculations can be simplified by employing the collective multipole expansion of the free 
Green function over the basis set of Wigner D-functions. The coefficients of that expansion 
are calculated in Sees. HV] IVl using two different approaches. In Sec. [IV] the expressions 
for multipole coefficients are derived by calculating the overlap integral of the free Green 
function with Wigner .D-function. The expansion of the free Green function over the set of 
six-dimensional hyperspherical harmonics is considered in Sec. [V] The derived results are 
discussed in Sec. [VI] where some concluding remarks are also given. Appendix [A] contains 
details of computations of the integral with the free Green function over the Euler angles. 
The summary of properties of the set of six- dimensional hyperspherical harmonics labeled 
by the particle's individual angular momenta quantum numbers is given in Appendix [B] 

The convention % = 1 is used throughout the text. Capital letters denote the six- 
dimensional vectors, e.g. R = (ri,r2), where 1*1, r 2 are three-dimensional vectors, R = 



|R| = ^r( + ri 



II. THE GENERAL FORMALISM 

We begin by writing the Schroedinger equation for the system of three particles with 

masses mi, m 2 , m 3 , 

fid 2 Id 2 \ 

(-^M-^*l + c, - £ )* (r " r2H0 - (1) 

Here, v x and r 2 are Jacobi vectors (see fig. d]) and the reduced masses /ii,/i 2 are defined by 
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FIG. 1: Jacobi vectors for the three-body system. CM23 is the CM of the particles m 2 and m^. 
It is convenient to replace r x — > J]T[yi, r 2 — > in eq. ([!]) which thereby becomes 



(A 6 -2p7- J E?l)*(r 1 ,r 2 )=0, A 6 = ^ + |1 (3) 

In this equation the potential energy depends on the masses of particles, but the kinetic 
energy operator does not. 

Below the situation is considered when the potential energy of the three-body system is 
much smaller than the total (non-negative) energy E. Introducing the wave number Q of 
the continuum state by E = Q 2 /2, we write 

(A 6 -2f/ + Q 2 )^( ri ,r 2 ) = 0, (4) 

We approximate the total wave function by the combination \1/ = + where \E f( - ^ is 
independent of U, and \E f< - 1 * ) scales as ^ — > c\I/ at U — ► cU . Substituting this decomposition 
into eq. (jlj) and comparing the coefficients which scale equally, we obtain 

(A 6 + g 2 )vI/(0)(r 1 ,r 2 ) = 0, 

(A 6 + Q 2 ) #W( ri ,r 2 ) = 2[/#(°). 

In the zeroth-order approximation the wave function is the product of two three-dimensional 
plane waves 

^ )(r 1 ,r 2 ) = e^ R ), (Q-R) = (q 1 -r 1 ) + (q 2 -r 2 ), (6) 

where qi i2 are the linear momenta vectors conjugated to the Jacobi vectors ri j2 . Conse- 
quently, the second equation in §5§ becomes 

(A 6 + Q 2 )¥ 1 \r 1 ,r 2 ) = 2Ue i ^ n \ (7) 

In order to solve this equation we first consider its homogeneous form 

(A 6 + Q 2 )$(R) = 0. (8) 



We are interested in a solution depending only on the hyperradius R. Therefore, we have 
to omit the angular part of the Laplacian Ag, which yields 



(9) 



Introducing the function F = $/i? 2 , this equation transforms to 



(10) 



The general solution of this equation is the combination of two Hankel functions (QR) 
and H¥\QR) 
spherical wave 
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. We note that H2 (QR) at large R has the asymptote of the outgoing 



(11) 



Thus, the solution of (El) can be chosen to be 



H^(QR) 
R 2 



$(R) 



(12) 



This function satisfies eq. (jHJ) everywhere except R = 0. Below we shall prove that 



G(R) = G(R) = -i 



Q 

AtiR 



H?(QR). 



(13) 



is the Green function of the equation (jSj), that is the solution of the inhomogeneous equation 



(A^ + Q 2 ) G(R) = 5(R). 



(14) 



Indeed, at R > this function satisfies the homogeneous equation (jSJ). At the limit R — > 
there is a divergency, 



G(R) 



Q 



(15) 



From the properties of Dirac ^-function it follows that integral of (1141) taken over the ar- 
bitrary region which includes the point R = must be equal to unity. We calculate the 
six-dimensional integral of (Till) taken over the infinitesimal sphere V e with its center in 
R = 0, 



(A 6 + Q 2 ) G(R)d 6 R = dS- V 6 G{R) + vr 3 Q 2 / G(R)R 5 dR. (16) 
\; JSe R=e Jo 



Here, the Gauss theorem has been used in order to transform the volume integral into the 
integral over the surface of 6-sphere S e with the radius e. Noting the expression for the 
vector surface element dS = (R/-R) e 5 dS, the integral (TTB|) at e — > evaluates to 

/ (A 6 + Q 2 ) G(R)d 6 R m [ dS - % f RdR = 1 - ^ -> 1. (17) 
Jv ( ^ Js e 4 Jo 8 

As is seen, G(R) defined by eq. (jT3l) does satisfy the equation (fl4l) . 

Having Green function (lT3i) . the solution of the equation ([7]) can be immediately written 

as 

^ (1) (R) = 2 J e l(Q - R,) C/(R') G(|R - R|) c/ 6 i?'. (18) 

Thus, the continuum wave function ^(R) of three particles within the perturbative approach 
has the form 

*(R) = *q(R) = e i «*)-i^ f e<^U(R')^^--^d e R'. (19) 

III. THE CALCULATION OF THE MATRIX ELEMENTS WITH THE CORREC- 
TION TERM 

In many physical applications it is necessary to calculate the matrix element 

M m = (tf^lOlttq), (20) 

where O is the transition operator and ^ is the wave function which corresponds to the 
state having the total angular momentum numbers jfi. Without loss of generality, we can 
assume that O is a scalar operator. This situation realizes, for example, in the case of the 
three-atomic molecular predissociation where O represents the operator of non-adiabatic 
couplings [5j. If the transition operator is a tensor Oi m then the product ^/j^Oi m can be 
decomposed into irreducible parts [9] so that the matrix element will be decomposed into 
the sum of matrix elements of the kind (1201) . 

We note that ( 120|) is the six- dimensional integral, 

M m = j ^(R)OV Q (R)d 6 R. (21) 

Among the six variables three can be chosen to be the collective angles determining the 
orientation of the whole system in space. The remaining three "shape" variables £ determine 
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the internal dynamics of the system and they can be chosen to be £ = r%, r^, cos#, where 9 is 
the angle between ri and r 2 . Thus, the transition operator O is the function of £, O = 0(£). 
Substituting the decomposition (TIT)]) into (12~T|) one arrives at the expression 

M, M = M^+M«, (22) 

where the zero- and first-order matrix elements are defined by 

AfjJ? = J ** M (R) Oe i(Q R,) d 6 #, 

M}2=-i^J **„( R) | e«**0 17(0 « ^ 



Here, by writing U{£') we make the assumption that the potential energy depends only on 
the three shape variables. This is true if the three-body system in not a subject of external 
forces. 

The calculation of the matrix element M„-„ has been considered in 71. Therefore, below 
we concentrate on the problem of calculation of the first-order (with respect to the potential) 



matrix element M]„ . 

As is seen, is determined by the 12-dimensional integral. However, six of twelve 

variables are the collective angles describing the rotation of the whole system from the 
coordinate frame defined by Jacobi vectors (ri, r 2 ) to the frame defined by (r' 1; r' 2 ) and from 
the frame (ri, r 2 ) to (qi, q 2 ). Noting that the potential energy does not depend on collective 
angles, one can try to integrate them out analytically. In order to do so one has to expand 
the integrands over the basis of angular functions depending on the collective angles. It is 
convenient to choose as the angular basis the set of Wigner D- functions 9j. 

The initial state wave function ^/j M (£) can be decomposed into the combination of (2j + 1) 
"internal" wave functions ipju(0-> depending on three shape variables £, 

* iM (R)= j^^ v {t)Di^), (24) 

u=-j 

where Q denotes three collective Euler angles describing the rotation from the body-fixed 
frame (BF) defined by the Jacobi vectors R = (ri, r 2 ) to the laboratory frame (LF) defined 
by the momenta vectors Q = (qi, Q 2 ). 

At this stage we have to evaluate the angular integral 

/ Q (R) = / eW) "g |R R T 2 R|) (25) 



where Q' denotes three Euler angles describing the rotation from BF' defined by the pair of 
Jacobi vectors R' = (r' 1; r' 2 ) to LF. 

The calculation of /q(R) can be performed by taking the multipole expansion of both 
integrand functions, 

oo j 



e *(Q.R') = £ Fftfa&D^W, 

(26) 



R'-R 2 ^ 8tt 2 

3=0 H,v=-j 

where Q" describes the rotation from BF defined by (i"i,r 2 ) to BF' defined by (r^rQ. 
In the upper equation £ g denotes three "shape" variables in the momentum space, e.g. 
£q = qx, q2, cosx, where x is the angle between qi and q 2 . 

Note that the rotation fl" can be presented as a product of rotations 

n" = Q x ft'' 1 . (27) 

In terms of Wigner D-functions this equation reads 

D U Q ") = E D iA Q ) D i'A^) = E D lA n ) { D -uM Q ')T- (28) 

v'=-j v'=-j 

Substituting this equation and eqs. ( 1261) into eq. ( |25i) . we arrive at the identity 
2/ + 1 



x | ^(^')[^ ) -,»(^)]*^ / - (29) 
The integration gives Sj 7 j'6/j,- u '8 v - v ti 8n 2 / (2j + 1) so that 

oo 

j q(R) = E E F< il^ o ^(£; o ^(«)- (so) 

j=0 ^jyu'.f 



Inserting this equation together with eq. (jMj) for the initial wave function into eq. f[23|) for 
the matrix element M-*', one can integrate out the angles f2. The result is 



i E / vkcewe) E / ^'(^0^(e ff ,o^(0^^- (31) 



m 2j + 

v=-3 It' =-3 



As is seen, the 12-dimensional integral fl23|) is reduced to the 6-dimensional integral. 
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In the next sections we consider the problem of the calculation of the collective multipole 
coefficients G^l(£; £') introduced in eq. fl26l) . As a preliminary, we note that these coefficients 
obey the symmetry relation 

[G^or = (-iy-"G^i_ u (t;o (32) 

which follows from the symmetry property of D-functions (Vt) = {—1)^ U [D 1 _ I1 _ U {VL)\*. 
Thus, in practice it is enough to calculate only the coefficients with /i> v. 

For the sake of completeness we present also the expression for the volume element of the 
shape space for two most often used choices of the set of shape variables £ 

r r°° f°° r f°° r /2 sin 2 2a r 

d^= / r\dr x / r\dr 2 / sm6d9= / R 5 dR / — — da / sm6d6, (33) 
J Jo Jo Jo Jo Jo ^ Jo 

where the second equation corresponds to the hyperspherical set £ = R,a, 9, where R = 

^Jr\ + r 2 ; is the hyperradius and a = arctan(r 2 /r!) is the hyperangle, so that 

r\ — R cos a, r 2 = R sin a, < a < tt/2. (34) 



IV. MULTIPOLE COEFFICIENTS FOR THE ANGLE-BISECTOR GAUGE 

Using the orthogonality of D-functions the expression for the multipole coefficients G^} v 
can be written as the three-dimensional integral, 

GgUfcO = / ^0^^\DU^Td 3 n''. (35) 

Here, apart of D-functions, Euler angles Q" enter the integrand through the term |R — R'| 
which is 

|R-R'| = ^R 2 + R' 2 -2p, p=(r 1 .r / 1 ) + (r 2 -r / a ), (36) 

where only the parameter p depends on Q". 

At this stage, one has to specify explicitly how the axes of the BF and BF' are connected 
to the Jacobi vectors (ri,r 2 ) and (r^r^). This procedure is not unique and, therefore, it 
can be denoted as the "gauge convention" [1 Of ] . 

Below we use angle-bisector gauge for both BF and BF', see fig. EJ Since p is scalar, it 
does not matter in which frame we calculate it. Let it be BF'. In this frame the components 
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FIG. 2: The angle-bisector gauge: the x-axis is directed along the bisector of the angle between 
the unit vectors fi and ?2- The z-axis is directed along the vector product [fj x £2]. BF' is defined 
analogously. 

of R' have the form 

r;=r;[cos(£72),-sin(#72),0], 
r' 2 =r 2 [cos(072), 8*11(072), 0], 

The components of R in BF' can be calculated using the equation 

(fn)<= J] {r n )kA ki {tt"), n = l,2. (38) 

n 

where *4fcj(fi") denotes the components of the rotation matrix [9J and (r n )fc is fc-th component 
of r n in BF where it is given by eq. (I37p with dashes removed. From eq. (1381) follows the 

identity p = Y,i,k(( r 'i)i + (r' 2 )i (r 2 )/cAi)- 

Using explicit form of the rotation matrix A, after some algebraic transformations de- 
scribed in detail in [71], the parameter p can be written as 

p — di sin(a + 7 + <5i) — a 2 sin(a — 7 + 5 2 ). (39) 

where the parameters are defined by 



pY ( . p x 2 



a x = px [ cos - l , 01 = p 2 I sm - 



(40) 



Pi,2 = ( r i r i) 2 + ( r 2r 2 ) 2 + 2r 1 r 2 r[r , 2 cos(0 =f 0'). 

where px corresponds to (0 — 0'). The additional phases Sx,2 are defined by 

rxr\ + r 2 r 2 T 0' 

sindi 2 = cos — - — , 

Pi,2 2 

r 2 r 2 — r[ri 9 9' 

cosdi 2 = sin — - — , 

Pi,2 2 



(41) 
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where the upper sign corresponds to 8\. 

In order to calculate the integral in fl35l) we employ the series expansion (eq. (7.15.18) of 
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fill), 

(Q\R> - R|) _ A (Qp) n Hi% (QVR 2 + R> 2 ) 



n! (m -L R/2W2+1 • ^ > 



IR'-RI 2 ^ n! (R 2 + R' 2 ) n / 2 + 1 

i i n=0 v ' 

Thus, the expression for the collective multipoles can be written as 



(P 2 + #' 2 W^ 

71=0 V V 

where the functions P„ ^ are defined by the integral 



(43) 



P "^ ) = ^! J P n l D U Q ")T d 3 n" . (44) 

This integral is calculated in Appendix |A] (see eq. (1A10I) ). In Appendix |A] it is also demon- 
strated that the functions Pn U ' are non-zero only for indices /i, v and n having the same 
parity. Hence, the multipoles & v vanish if p and v have different parity. 

Below we present explicit expressions for the parameter v for the most simple cases 
j = 0, 1. For 5-state the parameter pi° 00 ^ is given by eq. ( 1A12I) which leads to 

n=0 • 7 



where w = (p 2 + p|)/(2pip 2 ) and P n (w) denotes the Legendre polynomial. For P-state 
(j = 1) there are three multipoles: Gj , = — ((?!]_ -J*, ( j| = (G^ -J*. These 



coefficients correspond to the states with different spatial parity ll| . Namely, Gq is the 
pseudotensor parity multipole (state p even ^ and G\ ±l are polar tensor multipoles (state 
poddy p resen ^ fare only the expression for Gl , 

G ,oiSA )- R2 + R/2 { pip2 )2^ {2n + 2)] \ R 2 + R ,2) ** n w 

x H^ +2 (QVR 2 + R' 2 ) , (46) 

where P^iw) = dP n {w)/dw. 
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V. THE HYPERSPHERICAL EXPANSION OF THE MULTIPOLE COEFFI- 
CIENTS 



The hyperspherical form of the multipole coefficients can be derived based on the following 
expansion of the three-body free Green function 

^(QIR'-RI) " 

(ORR') 2 

n=0 



IR' - R| 



A 

(QRBjf X> + 2 ) ^ ■ J n + 2{QR) H%{QR!), R < R' , (47) 



where R • R' = p/(RR') and p is defined in ( 1361) ; C^(R • R') is n-th order Gegenbauer 
polynomial [8J]. If R > R' then the substitution R <->• R' must be made in ()4Tj) . 

According the definition (|35|) . the expression for the multipole coefficients can be written 

as 

1 oo 

= p^E^^^SM. r<#> ( 48 ) 

where J n+2 (QR) is Bessel function and the functions /ii"""^ are defined by the integral 

hgr) = 4 ( n + 2 ) y C^(R ■ R') [^(n*)]* (49) 
In order to calculate this integral we note that the Gegenbauer polynomial in (H7j) is pro- 



portional to the scalar product of six-dimensional hyperspherical harmonics 



it 



TV 2 



Cn{R ' R,) = 2FT2) {Yn{R) ' Fn(R,)) " (50) 

It is convenient to define harmonics Y n to be the eigenfunctions of the angular momentum 
operators corresponding to the Jacobi vectors ri and r 2 . In this approach, six- dimensional 
spherical harmonics are proportional to the product of conventional three-dimensional har- 
monics depending on the spherical angles of unit vectors fi and r 2 , 

Y nqlmm ,(R) = C^- 2 ^ l \a) Y lm (h) F (n _ 29 _ 0m ,(r 2 ), (51) 

where the functions Cq 2q ~ l,l \a) depending on the hyperangle a are defined by eq. (IB2I) of 
Appendix [B] 

The expression ( 1B1I) for the scalar product of six-dimensional harmonics contains the 
product of two Legendre polynomials depending on (fi ■ r[) and (r 2 ■ f' 2 ). Thus, in order to 
calculate the integral (j4"9j) it is necessary to decompose the product of two Legendre polyno- 
mials into a combination of angular functions corresponding to the states with well-defined 
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values of the total angular momentum j. This can be done by expressing the Legendre poly- 
nomials via three-dimensional spherical harmonics. Next, the product of spherical harmonics 
can be re-written in terms of bipolar harmonics, 

i+v 

• ? x ) JV(* 3 • f 2 ) = J2 (-^'-'(cf^hM-Cf'ii^)), (52) 
i=|J-i'l 

where the bipolar harmonics C"'(fi, r 2 ) are tensor product of two spherical harmonics 9J, 

Cj^(f 1 ,f 2 ) = ^Iml'm'^lrni^l) C i / m /(f 2 ). (53) 
m,m' 

Here, C\ m = a/47t/ (21 + 1) Y/ m are modified spherical harmonics and C^, m , are Clebsch- 
Gordan coefficients. The scalar product of bipolar harmonics does not depend on the choice 
of the coordinate frame in which it is calculated. Let this frame be BF', then 

V 

= £(-1)" ^ ^ i> f ' 2 )W [C!>(?i ; * 2 )W ^OT, (54) 

where the subscript BF (BF') denotes the coordinate frame in which the components must 
be calculated. The second identity in fl54l follows from the tensor transformation rule under 
the rotation of the coordinate frame [9|. Now the integral ( )49l) can be easily calculated using 
the orthogonality of D-functions, 

4-7T 3 

h (JH = ^y-j+v — J^(2l + 1) (2[n -2q-l] + l) C^- 2q -^\a) Cf - 2<? ~ M) (a') 

3 + ff,i=0 

x [^ir^°(fi, f a )]Bi* [Cjj^ft, f 2 )] BF , (55) 



where the indices g, Z take non-negative integer values so that n — 2q — I > 0. To derive the 
explicit form of hn it is necessary to specify the gauge convention. The only exception is 
the case of j = (S-wave part of the expansion flSD ), when the coefficient hn°°^ is scalar 



hT = 4vr 3 J> + 1) C#3(a) cS(a') P,(cos0) P(cos0')- (56) 



Note that the coefficients hn°°^ with odd n vanish, h^+i = 0. 
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For j > we use the angle-bisector gauge (see fig. [2] above). In that gauge, the components 
of bipolar harmonics in fl55l) have the form 



[C l £- 2q - l \h, U)]bf = C l /;- 2 "- l \n/2, -6/2, tt/2, 9/2) 

= J2 C lrn(n-2 q -l )(p - m) Cl m (n/2^ (57) 

m 

Note that the explicit expression for the spherical harmonics Ci m (ir/2, 9) is rather simple 
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CL[^,~) = r } V ( ' +m)!!( " m)!! ' / + m " eVen (58) 

- - 1 1 0, l + m = odd. 

From this identity and the symmetry properties of Clebsch-Gordan coefficients it follows 
that are non-zero only for n, [i and v having the same parity. 

For ?' = 1 and u = v = it is convenient to use the reduction formula for bipolar 

n 

harmonics |13| . 

fi(Mi)1 - vw/ w - (59) 

where P^(9) is the associated Legendre polynomial, Pl(9) = — sin 9 P( (cos 9) . As a result, 
the coefficient ( 1551) assumes the form 

hT = 4vr 3 it <&W Pl(9) Pt{6>). (60) 

Coefficients with odd index n vanish, /4n+i = 0- 



VI. CONCLUSION 



In the presented paper the perturbation theory was applied to the calculation of the wave 
function of the three-body system in which the potential U of the inter-particle interaction 
is small comparing to the kinetic energy. 

The expression for the wave function in the zeroth- and first-order approximation with 
respect to U is given by eq (1191) . The problem of the calculation of the matrix elements be- 
tween the continuum state wave function and the wave function having well-defined angular 
momentum quantum numbers has been considered in Sec. IIHI The calculation procedure 
leads to the appearance of twelve- dimensional integrals (see eq. (|23l) ) which, in the most 
important particular case of an isolated system, can be reduced to six-dimensional integrals, 
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see eq. ( )3T|) . The reduction was achieved by employing the technique of collective multipole 
expansions in terms of Wigner D-functions (see eqs. (1241) and (1261) ) describing the rotation 
from the body-fixed frame (BF) to the laboratory-fixed frame (LF). 

The collective multipole coefficients (CMC) of the expansion of the Green function of 
the system of three free particles were calculated in Sec. [IV] and |Vj The two different ap- 
proaches were used for the calculation of CMC. The first one (Sec. IIVP is the straightforward 
computation of the integral defining CMC using BF and LF corresponding to the so-called 
angle-bisector gauge (see fig. [2]). The expression for CMC in that approach is given by 
eq. 

The second (hyperspherical) approach was based on the expansion of the three-body free 
Green function over the basis set of six-dimensional h ype rspherical harmonics. Note that 
those were not conventional harmonics as defined in jg, 12] but the set of harmonics labeled 
by the individual angular momenta quantum numbers of particles [lij]. The expression for 
CMC was then extracted from the series by using some specific gauge convention (eqs. ( l48j ). 
(ESD of Sec. IV]). 

The advantage of the straightforward approach is that the resulting series representation 
of CMC is simpler compared to that of the hyperspherical approach. Also, it remains 
valid in the whole shape space while the hyperspherical series diverge at the configuration 
when R = R'. Nevertheless, the hyperspherical representation can be more favorable if the 
potential energy can be factorized in terms of hyperspherical coordinates. In this case the 
many-dimensional integrals of the perturbation theory reduce to the one- dimensional form. 

Both representations of CMC are series of Hankel functions Hn ' depending on hyperra- 
dial variables multiplied with weight functions depending on the hyperangles. These series 
converge rather well. For example, a.t Q = 10, R = 1, R = 2, a = a' = 9 = 9' = 45° the 
number of terms in series which are needed to achieve the accuracy of 10~ 6 is about thirty. 
The convergence rate decreases as Q and/or i?'s increase. The same behavior is observed as 
the configuration triangles (see fig. d]) become narrower at the same Q and i?'s. 

In general case of j > 0, CMC are functions of six variables, i.e. two times the dimen- 
sionality of the shape space of three particles. Surprisingly, from eq. f]4"5]) it follows that 
CMC corresponding to an S-state (with j = 0) depends only on three variables. (Note that 
this property of S'-state CMC is not seen from 
property of S'-state CMC was already noticed in 



56]) of the hyperspherical approach.) Such 
7| where the collective multipole expansion 
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of the product of three-dimensional plane waves has been considered. This fact still needs 
some physical explanation. 

Currently, the work on the application of the developed theory to the problem of the 
three-particle fragmentation process is in progress. 
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APPENDIX A: THE INTEGRATION OF p n OVER EULER ANGLES 



In order to calculate the integral (14411 we employ the expression for Wigner D-functions 
in terms of Euler angles, Dj^ = e~ t<ya+ ^ u ((3). First, we calculate the integral over the 
angles a and 7, 



/*27r r2ir 

l n ^= \ \ e^ Q+ ^ ) (a 1 sin(« + 7 + 5 1 )-a 2 sin(a-7 + 5 2 )) n rf«rf 7 , (Al) 
Jo Jo 

This integral can be calculated using the exponential representation of sine functions. It 
turns out that the integral is non-zero only for indices n, fi, v having the same parity. 
Omitting details of routine transformations, we write 

I n—u—2k 2k+v 
J -A^( Ilf r ipil p -gi2 n! V a l a 2 (A2) 

nilv ~ [ } ^ k\ (k + u)\ - k)\ - ky: [ ] 

where the sum runs over all non-negative values of k at which factorials remain finite. Note 
that factorials and indices £l, v in eq. (1A2I) are always integer numbers (as was mentioned 
above) , 

u + v 11 — V , . . 

t =!L 2-> 9 = ^~- ( A3 ) 

Using the explicit form of the function c?^(/3) we can integrate (1A2I) over f3 term by term. 

The partial integrals are 



a 







r r~ 2k 4 h+P dlM sin(3 d(3 = 2v/(j + (j - fj)\ (j + v)\ (j - v)\ 

y y n n-p-2k n 2k+v (~iy-^ ( q + n + v-2k)\(j-q-v + 2k)\ 

Y q\U-li-qW-u-q)\(u + ii + q)\(j + n + l)V 1 J 
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Substituting eq. (1A4|) into eq. (1A2I) we can re-write the summation over k in (1A2j) as 



p 2 y fc (g + n + i/-2fc)!(j-g-i/ + 2fc)! ^ 



Pl 



2k 



(2k + /i-u- 1)!! (n-fi-2k- 1)H (g + n + i/ - 2Jfe)! (J-q-u + 2k)\ 
X fc! (2±ii _ (n-//-2fc)!(/x-i/ + 2fc)! ' ( ' 

The last multiplicand on rhs on this equation can be presented in a compact differential 
form 

(q + n + v - 2k)\ (j-q-v + 2k)\ 



t =r (A6) 



{n - fi - 2k)\ (p - u + 2k)\ 
Substituting this identity into (IA5I) we put derivatives outside the summation. Thereby, it 
evaluates in closed form 



^(...) = (_l)M+-+? 2 ^ 



v(^-i/-l)!!(n-^-l)!! 



gj-H-q j.n+j+1 Qfi+v+q j.n-n-1 



' Pi 



1 + /i — i/ n + z/ 1 + \i — n (tp 2 ) 2 
~~2 ' 2~ ; ~ 

where 2 i*i denotes Gauss hypergeometric function fis| . 

The final step is to rearrange the exterior summation over q 

1 



, (A7) 



i=l 



E 



?! (j - P - ?)! (i - f - ?)! + P + ?)! 



x 



E 



(n + j + 1)! 



Qj+v-m 



Qj-V-Q fn+j+l Qq+fi+u _ ^n+j+1- 



1 



m! (n + j + 1 — m)! ^— ' q\ (j — u — q) \ {y + fi + q) \ (j — // — g — m) 



, (A8) 



Here, the decouple sum over q evaluates to a closed form, 
Vr , (2j-m)! 



(j + A*)! (j ~ f )! (j - P - ™) ! U + v ~ m )- ' 



(A9) 



Collecting above eqs. (1A2I) . (1A4I) . ( 1A7I) - (1A9I) we finally write the expression for the 
integral (EHl) . 



x 



(2j-m)! 



(j + /i)! (j — i/)! ^— ' m! (j — /j, — m) \ (j + is — m) \ (n + j + 1 — m)\ 

1 + /j, — v n + v 1 + fj, — n (tp 2 ) 2 



x ^' +y - m ^- n - 1 2 F 1 



' 2 
Pi 



(A10) 



t=i 
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Note that the derivative can be written as a combination of Gauss hypergeometric functions, 

1 + \i — v n + v 1 + fi — n {tp2) 2 



gj+u-m t fi-n-l a _p i 



Pi 



t=i 



a=(j+v— m), 
(j+v— m) — 1,... 



(2a — j — v + m) \ (j + v — m — a)\ 



xl lz£±it) , Fl (l±>pi,-^,l±Jf±- a ,ti h ( Aii) 



where ( - ^ +n ) a denotes Pochhammer symbol [1 



Using eq. ( 1A10I) . one can easily calculate the coefficient P^^ corresponding to an S-state, 



P, 



(000) 



.2 „2n 



(2n-l)!! 



2n 



(A12) 



(2n)!! (2n + 1)! 

Here, the hypergeometric function reduces to the Legendre polynomial which leads to 
eq. ( H5i) . For P e -state with pseudotensor parity the non-zero parameter is 

2 (PlP2) n 



P, 



(100) 



2n 



-16tT 



(2n + 2)! 



2tpip 2 



(A13) 



t=i 



Calculating the derivative in this identity we arrive at eq. f|46|) of the main text. Note that 
both above parameters with odd index vanish, i.e. P 2 n+i = An+i = 0- 



APPENDIX B: THE SIX-DIMENSIONAL HYPERSPHERICAL HARMONICS 
AND THEIR PROPERTIES 



For the sake of completeness below we summarize the most important properties of the 



set of six-dimensional hyperspherical harmonics which are eigenf unctions of the 
angular momentum operators. The derivation of the results can be found, e.g. in 
The scalar product of two hyperspherical harmonics is defined by 



oarticle's 



(F n (R) • y n (R')) = E E E Y nqlmml (R') = £ LL±- 



47T 

1=0 9=0 mm' ql 

x (2[n - 2g - Z] + 1) Cf ~ 2 ^ l \a) Cf V) ■ r[) P N ^(v 2 ■ f 2 ), (Bl) 



where q m ax is the integer part of the ratio (n — l)/2 and the hyperangle a is defined in ([3 
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The functions Cg n 2q are defined as 



of Hoc) 



i 2 W'+3 (2q + 1 + 1' + 2) q\ (q + 1 + 1' + V 



7T (2q + 21 + 1)!! (2q + 21' + 1)!! 

x 2 9 (sina)'' (cos a)' Pf +1/2 ' m/2) (cos2a), (B2) 

where p^ a,b ^ is Jacobi polynomial 8]. Above functions are orthonormal, 

^CfHa)C^\a)^da = 5 M , (B3) 



The hyperspherical harmonics defined by eq. (1511 are orthonormal, 

j ' Y nq [ mm /(R) Y nqlmm , (R) d n = d n ^ ni d qm di^ x d m ^ mi d m t !m ^, (B4) 

where <i 5 r2 is the surface elements of the six-dimensional hypersphere, which can be written 

as 

<i 5 fi = — - — sin Q\ sin ^2 da d0\ d02 d<\>\ d<p2-, (B5) 

where 8\ y 2, 4>i,2 are three-dimensional spherical angles of the unit vectors f 1^. Note that the 
total solid angle in the six-dimensional space is 

d 5 Q = vr 3 , (B6) 

s 6 

where 5*6 is the six- dimensional hypersphere of the unit radius. 
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